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In recent years, there has been increasing interest in the specific heat C of insulators and semi- 
conductors because of the availability of samples with different isotopic masses and the possibility 
of performing ab initio calculations of its temperature dependence C{T) using as a starting point 
the electronic band structure. Most of the crystals investigated are elemental (e.g., germanium) 
or binary (e.g., gallium nitride) semiconductors. The initial electronic calculations were performed 
in the local density approximation and did not include spin-orbit interaction. Agreement between 
experimental and calculated results was usually found to be good, except for crystals containing 
heavy atoms (e.g., PbS) for which discrepancies of the order of 20% existed at the low temperature 
maximum found for C/T^. It has been conjectured that this discrepancies result from the neglect 
of spin-orbit interaction which is large for heavy atoms (Ao ~1.3eV for the p valence electrons 
of atomic lead). Here we discuss measurements and ab initio calculations of C(T) for crystalline 
bismuth (Ao ~1.7 eV), strictly speaking a semimetal but in the temperature region accessible to 
us (T > 2K) acting as a semiconductor. We extend experimental data available in the literature 
and notice that the ab initio calculations without spin-orbit interaction exhibit a maximum at 
~8K, about 20% lower than the measured one. Inclusion of spin-orbit interaction decreases the 
discrepancy markedly: The maximum of C{T) is now only 7% larger than the measured one. Exact 
agreement is obtained if the spin-orbit hamiltonian is reduced by a factor of ~0.8. 



I. INTRODUCTION 

In the past few years a number of investigations on the 
dependence of the specific heat of semiconductors and in- 
sulators on temperature and isotopic masses have been 
carried out. These works involved careful low ternpera- 
ture experiments (for elemental crystals see Refs. [l| and 
H, for binaries see Refs. H and references therein) and 
elaborate ab initio calculations based on electronic band 
structures computed in the framework of the local den- 
sity approximation (LDA). Recent work involving binary 
lead chalcogenides with different isotopic compositions 
(PbS, (Ref. lil), PbSe (Ref. ^) yields a low temperature 
maximum in the quantity C{T)/T^ (at ~12K)ii^ about 
25% lower than the measured one. Correspondingly, the 
calculated phonon dispersion relations are on the average 
10% higher than those determined with inelastic neutron 
scattering (INS)^.The higher phonon frequencies qualita- 
tively explain the lower calculated specific heats. 

When trying to ascertain whether the harder phonon 



frequencies (and lower specific heats) calculated for PbS 
and PbSe were due to the lack of spin-orbit (s-o) cou- 
pling in the ah initio electronic structure some difficul- 
ties arose. These involved not only the extended compu- 
tational time but also divergences in the dispersion rela- 
tions of the optical phonons for g— >0. The latter may be 
related to the strongly ionic, near ferroelectric character 
of these compounds. 

To our knowledge the computation of the q^O di- 
vergences when the spin-orbit coupling is present has 
not yet been implemented in existing first-principle soft- 
ware. By contrast, for metals, there is no such diver- 
gence, which makes possible the direct computation of 
phonon band structure and thermodynamical proper- 
ties based on Density-Functional Perturbation Theory 
(Ref. 7,8,10), including spin-orbit, as implemented in the 
ABINIT software®. Thus, we performed measurements 
and calculations for crystalline bismuth, which is free of 
the ionic divergences present in the lead chalchogenides. 

In this Letter we demonstrate that inclusion of s-o cou- 
pling, in fact, considerably reduces the discrepancies be- 
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tween experimental heat capacity data and ab initio re- 
sults. In addition, we also discuss the dependence of the 
rhombohedral lattice parameter ao on the magnitude of 
s-o coupling. 

Bismuth is a semimetal closely related to the lead 
chalcogenides: It has 10 valence electrons per primitive 
cell and a rhombohedral structure which can be derived 
from that of PbS by making both atoms equal and ap- 
plying a Peierls-like distortion to the PbS cube, involving 
an elongation of one of its [111] axes^^. The distortion 
converts the simple cubic structure of bismuth, with one 
atom per primitive cell, into a rhombohedral one with two 
atoms per primitive cell (two sublattices). This structure 
is characterized by three parameters: the bond length 
ao, the rhombohedral angle and a shift between the two 
sublattices^ii^. These parameters were determined by 
energy minimization using the ABINIT code. The re- 
sults obtained for these parameters in Refs. [T^ and [l^ 
with and without s-o coupling differ by less than 1%. 

The phonon dispersion relations were calculated with 
ABINIT, in Ref. [ij with s-o coupling whereas in Ref. 
[13I calculations with and without s-o coupling were per- 
formed. With s-o coupling, excellent agreement with 
experimental results was obtained whereas without it, 
discrepancies of the order of 10% were found, the cal- 
culated bands lying higher in frequency than the mea- 
sured ones. Similar results were found for PbS (Ref. i) 
and PbSe (Ref. i) without s-o coupling. In view of 
these results we proceeded to calculate the specific heat 
C{T) of bismuth. Measurements which were performed 
on high purity (99.9999%, Preussag Pure Metals) single 
crystals in order to complete the available experimental 
data (performed for polycrystalline samples 
Recent data collected by various earlier authors and our 
own data (the experimental technique is described in 
Refsdlilllj) are compiled in Fig. [TJ together with the 
results of our ab initio calculations performed with and 
without spin-orbit coupling. 



II. RESULTS AND DISCUSSION 

The hitherto available experimental points were mea- 
sured on polycrystalline samples. They are rather widely 
spaced in temperature, with the exception of Keesom's 
below 4 K (RefLl6''20). We have therefore performed 
measurements of Cp{T) on single crystals for T between 
1.8 and 100 K with the measuring temperatures spaced 
by -O.IK up to 50K and 0.5K steps above 50K. The 
maximum oi C{T)/T^ takes place at 7.5 K. According to 
Ref. i it should be found at - Tta/Q, where Tta -^42 K 
is an Einstein oscillator frequency which can be read off 
the phonon density of states^. The temperature of the 
maximum in Fig. [l]is found to be 7.5 K, fairly close to 
Tta/6 = 7 K. 

Our measurements place the maximum value of C/T^ 
at 2.320 ± 0.03 mJ/mol K"*. The value calculated with- 
out s-o splitting lies at 1.940 mJ/mol K** whereas with 
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FIG. 1: (Color online) Heat capacity of a Bi single crys- 
tal, purity 99.9999%. (red filled) circles as measured in this 
work compared with literature data obtained on polycrys- 
talline samples, (green) ▲ (Ref. [H); (blue) ♦ (Ref. 'l?'); □ 
(Ref. [H); (black) ■ (Ref. M)- (black) solid line: ABINIT 
results with spin-orbit coupling included; (blue) dashed line: 
ABINIT calculation without spin-orbit coupling. 



s-o coupling one finds 2.500 mJ/mol K^. Hence, the cal- 
culation without s-o coupling lies 27% below the experi- 
mental data, that with s-o coupling only 7% above. Not 
only is the difference between measured and calculated 
values of C /T^ three times smaller when s-o interaction is 
taken into account, but the agreement with s-o coupling 
above 12 K lies within the experimental error, whereas 
without s-o coupling a considerable difference is found. 
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FIG. 2: (Color online) Calculated heat capacities of Bi with 
varying magnitude of the s-o coupling as indicated in the in- 
set. Note that the maximum shifts to lower temperatures 
(inclined dashed line) with increasing s-o coupling, as corre- 
sponds to decreasing phonon frequencies (see text). 



3 



This improvement illustrates the importance of s-o in- 
teraction for the calculations of thermodynamic proper- 
ties starting from ab initio electronic band structures for 
systems containing heavy atoms and confirms our con- 
jecture that similar discrepancies between calculated and 
measured values of C/T^ found for PbS and PbSe^i^ are 
also due to the lack of s-o coupling in the electronic 
structure calculations. These results suggest perform- 
ing calculations and measurements for antimony (also a 
semimetal with the same crystal structure as Bi) which 
has a considerably smaller atomic s-o coupling than bis- 
muth (0.84 vs. 1.7eV). 

While the Sb work is in progress we have pursued yet 
another avenue: we have multiplied the spin-orbit cou- 
pling Hamiltonian by a factor < A < 1 and repeated 
the full ab initio calculations of C(T) for several values 
of A. The results obtained for bismuth are shown in Fig. 
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FIG. 3: (a) Maxima of the quantity CpiTmax) /Tmax vs. spin- 
orbit coupling parameter A. o our data, dashed line: fit with 
eq. [T]with parameters given in the text, (b) Energy minimized 
lattice parameter ao vs. spin-orbit coupling parameter A. o 
our data, dashed line: fit with eq. [l]with parameters given in 
the text. 



This figure clearly reveals the strongly supralinear de- 
pendence of the s-o effect on C /T^ which can bet fitted 
with the expression 



C/r3^co[l + c2A2(l + c4A2)] (1) 
with co=l. 942(2) mJ/molK'', C2=0.116(7) and 
C4=1.46(13). 

The fit is displayed in Fig[3I^a). For A=l the fourth or- 
der term surpasses the quadratic term. Exact agreement 
between the measured and the calculated values of C /T^ 
is obtained if one reduces the s-o coupling by the factor 
A=0.8. 

It is expected that other mechanical and thermody- 
namic properties of Bi should depend on the s-o coupling 
as well, i. e., on the value of A. The simplest of them 
is probably the rhombohedral lattice parameter oq. It 
is shown in FiglS^b) calculated for the same values of 
A as C{T) and fitted with an expression similar to eq. 
ID with ao(A = 0) = co=4.6529(7)A, 02=0.0069(9) and 
04=0.29(16). The values of the fit parameters q indicate 
that there is no simple relationship between the effect of 
s-o splitting on C{T) and the lattice parameter ao(A). 
The effect of spin-orbit interaction on C{T), as described 
by C2 and C4, is much larger than that on the lattice pa- 
rameter ao(A). For A=l, however, the calculated value 
for ao (4.694 A) is also much closer to the the experi- 
mental one (4.720 A) than that obtained for A=0 (4.652 

A). _ 

In conclusion, we have investigated the effect of s-o in- 
teraction on two thermodynamic properties, C(T) and oq 
of a solid consisting of the heavy element bismuth. Such 
effects, apparently rather substantial for a first-principles 
calculation of the physical properties, have not received 
much consideration in the literature, so far. For a dis- 
cussion of other similar phenomena see also Refs. [l^ and 
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